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ABSTRACT 

We investigate the effect of magnetic fields on the propagation dynamics 
and morphology of overdense, radiatively cooling, supermagnetosonic jets, 
with the help of fully three-dimensional smooth particle magnetohydrodynamic 
simulations. Evaluated for a set of parameters which are mainly suitable for 
protostellar jets (with density ratios between the jet and the ambient medium 
f] ^ 3 — 10, and ambient Mach number Ma ~ 24), these simulations are also 
compared with baseline non-magnetic and adiabatic calculations. Two initial 
magnetic field topologies (in ~ equipartition with the gas, f3 = pth/pB — 1) 
are considered: (i) a helical held and (ii) a longitudinal held, both of which 
permeate both the jet and the ambient medium. 

We hnd that, after amplihcation by compression and re-orientation in 
nonparallel shocks at the working surface, the magnetic held that is carried 
backward with the shocked gas into the cocoon improves the jet collimation 
relative to the purely hydrodynamic (HD) systems, but this ehect is larger in 
the presence of the helical held. In both magnetic conhgurations, low-amplitude, 
approximately equally spaced (A ~ 2 — 4i?j) internal shocks (which are absent 
in the HD systems) are produced by MHD Kelvin-Helmholtz rehection pinch 
modes. The longitudinal held geometry also excites non-axisymmetric helical 
modes which cause some beam wiggling. The strength and amount of these 
modes are, however, reduced (by ~ twice) in the presence of radiative cooling 
relative to the adiabatic cases. Besides, a large density ratio, rj, between the jet 
and the ambient medium also reduces, in general, the number of the internal 
shocks. As a consequence, the weakness of the induced internal shocks makes 
it doubtful that the magnetic pinches could produce by themselves the bright 
knots observed in the overdense, radiatively cooling protostellar jets. 



Magnetic fields may leave also important signatures on the head morphology 
of the radiative cooling jets. The amplihcation of the nonparallel components 
of the magnetic helds, particularly in the helical held geometry, reduces the 
postshock compressibility and increases the postshock cooling length. This 
may lead to stabilization of the cold shell of shocked material that develops 
at the head against both the Rayleigh-Taylor and global thermal instabilities. 
As a consequence, the clumps that develop by fragmentation of the shell 
in the HD jets tend to be depleted in the helical held geometry. The jet 
immersed in the longitudinal held, on the other hand, still retains the clumps 
although they have their densities decreased relative to the HD counterparts. 
As stressed before (Cerqueira, de Gouveia Dal Pino & Herant 1997), since the 
fragmented shell structure resembles the knotty pattern commonly observed in 
HH objects behind the bow shocks of protostellar jets, this result suggests that, 
as long as (equipartition) magnetic helds are present, they should probably be 
predominantly longitudinal at the head of these jets. 

Subject headings: ISM: jets and outhows - MHD - stars: pre-main- 

sequence - stars: mass loss 
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1. Introduction 

There is increasing evidence that protostellar jets are driven by circumstellar 
magnetized disks associated with pre-main sequence stars (e.g. Konigl & Ruden 1993; 
Shu et ah 1995). Efficient mass loss in supersonic, collimated magnetized outflows is the 
most likely mechanism by which protostars dissipate the angular momentum accumulated 
during the accretion of the surrounding material. While of fundamental importance in the 
production and initial collimation of the jets, magnetic fields have been generally neglected 
in most of the analytical and numerical modeling of the structure of protostellar jets since 
the inferred estimates of their strength {B rsj 10 ® — 10 ® G) suggested that they could 
be not dynamically important along the flow (e.g., Morse et al. 1993). However, recent 
observations of circularly polarized radio emission of the young stellar outflow of T Tauri 
S (Ray et al. 1997), for example, indicate the presence of a strong, ordered magnetic held 
in the flow, far away from the source, which has possibly been amplified by compression 
behind the shocks at the head of the outflow. After amplification and re-orientation behind 
the shocks, such magnetic fields may operate significant changes on the dynamics and 
collimation of a jet and also in its head structure as suggested by recent numerical studies 
(see below). 

Great effort has been concentrated in the analytical and numerical study of 
magnetized, adiabatic, light jets, aiming mostly the investigation of extragalactic jets (see, 
e.g. Birkinshaw 1997 for a review). Most of that work has focused on the study of the 
stability properties of the beam against hydromagnetic and Kelvin-Helmholtz instabilities. 
In the limit of zero-velocity difference between the jet and the surrounding medium, linear 
theory predicts that a jet magnetically confined by a toroidal held is unstable to the pinch 
and kink (or helical) hydromagnetic modes (e.g., Ghandrasekhar 1961) with the temporal 
growth rates of the pinching mode increasing with increasing density ratio, t], between the 
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jet and the ambient medinm. In the presence of a non-zero velocity discontinnity at the 
bonndary layer separating the two fluids, these pure hydromagnetic modes are modihed by 
the development of the Kelvin-Helmholtz (K-H) instability. A jet conhned by a toroidal 
magnetic held is unstable to the fundamental and rehection pinch and kink modes of the 
K-H instability (Cohn 1983; Fiedler & Jones 1984). The most unstable pinching mode in 
this case, at wavelengths A ~ 2 x iiRj (where Rj is the jet radius), has a destabilization 
length I oc {MjRj) (where Mj = vj/cj is the jet Mach number), which is similar to the 
pure hydrodynamical (HD) case for jets with r] < {Mj/2)^'^ (Cohn 1983). In the case of 
magnetized jets with a longitudinal magnetic held residing in an unmagnetized medium, all 
K-H modes become stable for sub-Alfvenic hows (except for a small region of slow rehection 
modes; e.g., Bodo et al. 1989; Hardee et al. 1992; Hardee, Clarke & Rosen 1997). In 
the super-Alfvenic regime, on the other hand, the growth rates of the instability are not 
very much diherent from those of pure HD hows. For example, a fundamental kink mode 
in the pure HD case can be identihed with an Alfven disturbance of long wavelength in 
the magnetohydrodynamic (MHD) case, and rehection modes at shorter wavelengths can 
be identihed with fast magnetosonic (ms) waves rehecting oh the jet boundaries. In this 
super-Alfvenic regime, if the jet is also supermagnetosonic {Mms,j = Vj/{v\ -f- c|)^/^ > 1, 
where va is the Alfven velocity and Cj is the jet sound speed), then it becomes more stable 
with increasing with the destabilization length varying approximately proportional 

to Mms,j {I oc MmsjRj) (Hardee et al. 1992). Also, strong toroidal helds of strength 
comparable to the poloidal held can lead to increased jet stability (Appl & Camenzind 
1992; Thiele & Camenzind 1997). 

The ehects of the K-H modes on the survival of the beam, however, cannot be predicted 
by the linear theory alone. As in pure HD hows, the fastest growing, shortest wavelength 
rehection modes are expected to saturate with the formation of weak oblique shocks, while 
the fundamental modes may not saturate and may cause large-scale distortions and even 
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disruption of the flow. Numerical simulations which can determine the end points of the 
operation of these instabilities conhrm these predictions. Hardee et ah (1992) and Hardee, 
Clarke & Rosen (1997), for example, assuming slab and cylindrical jets, respectively, with 
axial magnetic helds, have focused on the comparison of the scale-length of the structures 
generated during the nonlinear evolution and the wavelengths of maximum growth rate 
predicted by the linear theory. They hnd that the jet is not stabilized by nonlinear effects 
associated with increasing magnetic tension and disrupts near the resonant wavelength of 
the K-H kink mode. Malagoli, Bodo & Rosner (1996) and Min (1997a, b) have extended 
those investigations of the nonlinear development of the K-H instability by including 
diffusion and magnetic reconnection effects and found that even if the magnetic held 
intensity is not too large to completely suppress the K-H instability, it is still able to 
mediate turbulence decay and diffusion of energy and mass across the boundary layer 
between jet and surrounding medium. 

These studies of the stability properties of the beam against MHD K-H pinch and 
helical modes seem to provide potential mechanisms to explain the formation of structures 
such as knots and wiggles in adiabatic, supermagnetosonic light jets. Besides these stability 
analyses, some effort has also been spent in numerical studies of the general effects of 
B-£elds on the global evolution and morphology of adiabatic, light jets, still in the context 
of extragalactic jets, assuming both passive (Clarke, Norman & Burns 1989; Matthews 
& Scheuer 1990; Hardee & Norman 1990) and dynamical important toroidal fields (e.g., 
Clarke , Norman & Burns 1986; Lind et ah 1989; Kossl, Muller & Hillebrandt 1990). In the 
later case, the jet was found to be rapidly decelerated at the Mach disk with the shocked 
jet material being confined to a slender trans-Alfvenic plug instead of being deposited in 
the large cocoon observed in the purely HD light jets. 


Lately, these numerical MHD studies have been extended to heavy, adiabatic jets 
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(Todo et al. 1993; Hardee & Clarke 1995; Stone, Xu & Hardee 1997). Hardee & Clarke 
and Stone et al. have focused on simulations of jets with a poloidal held propagating into 
an unmagnetized medium, while Todo et al. (see also Thiele & Camenzind 1997) have 
considered a helical held conhguration extending also to the ambient medium. Comparing 
the development of supermagnetosonic, adiabatic, heavy jets with axial and toroidal 
magnetic helds propagating into unmagnetized ambient media, Hardee & Stone (1997) have 
found that the toroidal geometry suppresses the mixing and entrainment of ambient gas 
which is found to develop in the axial case as a consequence of the K-H instability. Also, 
they have found that the kink mode has a longer wavelength and smaller amplitude in the 
toroidal conhguration. 

Still almost unexplored, however, is the role played by B-helds on the propagation 
dynamics and morphology of radiatively cooling, heavy jets (a scenario which is believed 
to occur in protostellar jets). In the limit of zero magnetic held, numerical simulations of 
radiatively cooling, heavy jets [e.g., Blondin, Fryxell and Konigl 1990, hereafter BFK; de 
Gouveia Dal Pino & Benz 1993, hereafter GB93; Stone & Norman 1993a, 1993b, 1994; 
de Gouveia Dal Pino & Benz 1994, hereafter GB94; Ghernin et al. 1994; de Gouveia Dal 
Pino, Birkinshaw & Benz 1996, hereafter GBB96; de Gouveia Dal Pino & Birkinshaw 1996 
(GB96)] have shown that thermal energy losses by the jet system have important ehects 
on its dynamics. These studies have revealed that a cooling jet develops a dense, cold 
shell of shocked material at the head which is fragmented into clumps by Rayleigh-Taylor 
instability. Moreover, BFK and GB93 have found that the development of the K-H modes 
along the beam are inhibited by the presence of cooling. Recently, Hardee & Stone (1997) 
and Stone, Xu & Hardee (1997) (see also Massaglia et al. 1992) have examined the dynamics 
of K-H unstable cooling jets. Their linear analysis indicate that the growth of the K-H 
modes is very sensitive to the assumed form of the cooling function. In particular, if the 
cooling curve is a steep function of the temperature in the neighborhood of the equilibrium 



state, then the growth of K-H modes is reduced relative to the adiabatic jet - a result which 
is consistent with previous numerical simulations (BFK; GB93). With the inclusion of a 
longitudinal magnetic held in a supermagnetosonic jet, Hardee & Stone (1997) hnd that 
the magnetic held does not strongly impact the diherences in the K-H stability properties 
between an adiabatic and a cooling HD jet, provided that the magnetic pressure does not 
dominate. Besides, they hnd that the increase in the magnetic held strength makes the 
linear stability properties to become more like those of an adiabatic jet. The nonlinear 
analysis of the growth of the K-H modes in cooling jets in the 5 = 0 limit shows similar 
behaviour as in the pure adiabatic cases. The jet can disrupt near the resonant frequency 
of the fundamental kink K-H mode, and nonlinear, higher frequency, rehection waves tend 
to produce low-amplitude wiggles, and can result in strong shocks in the jet beam. 

In the present work, we attempt to extend these previous investigations by exploring 
the nonlinear ehects of magnetic helds (close to equipartition with the gas) on the global 
evolution and morphology of radiatively cooling, heavy jets and test the morphological 
signatures of two different magnetic held geometries (a longitudinal and a helical 
conhguration) on the dynamics of protostellar jets. A preliminary step in this direction was 
made in previous work (de Gouveia Dal Pino & Gerqueira 1996; Gerqueira, de Gouveia Dal 
Pino & Herant 1997, hereafter Paper I), where we have mainly focused on the effects of the 
magnetic helds in the jet head structure. With the help of three-dimensional (3-D) smooth 
particle magnetohydrodynamical (SPMHD) simulations, we have found that the presence 
of a helical magnetic held (in close equipartition with the gas) may suppress the formation 
of the clumpy structure that is found to develop at the head of HD jets by fragmentation 
of the cold shell of shocked jet material. A cooling jet immersed in a longitudinal magnetic 
held, on the other hand, tends to retain the clumpy morphology at its head. In the present 
work, we perform 3-D SPMHD simulations to address the details of the magnetic held 
ehects over the whole structure of the radiatively cooling, heavy jets, covering a more 
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extensive range of parameters, and compare with both non-magnetic and adiabatic systems. 

Contemporaneously with this work, Frank et al. (1998; see also Frank et al. 1997) have 
performed grid-based 2|-D simulations of magnetized, radiatively cooling jets assuming 
a toroidal {B^) magnetic held and concluded that toroidal helds are able to excite the 
development of strong pinches along the beam. Although in our simulations we have 
assumed somewhat different initial conditions (see below), where both investigations 
overlap, the results are qualitatively similar, except for the fact that in our simulations 
the MHD pinch modes are found to be not strong in radiatively cooling, heavy jets. This 
difference, however, is mainly due to differences in the assumed initial density ratios between 
the jet and the ambient medium, which are much larger in our analyses (see Paper I and 
discussion below). 

In §2, the numerical method and the initial conditions used are briefly described. In 
§3, we present the results of our simulations, and in §4, we address the conclusions and the 
possible implications of our results. 


2. Numerical Technique 

In our simulations, we consider the magnetohydrodynamics (MHD) conservation 
equations in the ideal approximation 


dp 

dt 


—pS/ ■ V 


(la) 


dv 

dt 


_Vp ^ X B) X B 

p Airp 


( 16 ) 
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§ = (ic) 

at p 

pjy 

_ = _B(V ■ v) + (B • V)v {Id) 

where the symbols have their usual meaning (i.e., p is the density; B is the magnetic held; 
u is the specihc internal energy, C is the radiative cooling rate, etc.). To close the system 
( 1 ), an ideal equation of state is assumed: 


P = (7 - 1)PM 


(le) 


with 7 = 5/3. 

The MHD equations above are solved (in Cartesian coordinates) using a modihed 
version of the fully three-dimensional smooth particle hydrodynamics (SPH) code originally 
developed by de Gouveia Dal Pino & Benz (1993) (see also GB94; Chernin et al. 1994; 
GB96; GBB96) for the investigation of the evolution of purely hydrodynamic (HD) jets. 

In the SPH formalism, equations (1) are described by (e.g., Benz 1990; Stellingwerf & 
Peterkin 1990; Monaghan 1992; Meglicki 1995): 



^ + ^ + n.. 




+ 


1 ^ 


Bj) X ViWij] X Bi 


(2a) 



Pi ^ 1 ^ 

^ H - ^j)'^iWij + 7 ; mjUij{Yi - Vj) ■ ViWij 

Pi j=l ^ j=l 


m 
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V i^Vij 


( 2 c) 


where mj is the mass of the particle j (located at the position r = rj); pi and pj are the 
density of the particles i and j, respectively; pi and pj are their pressure; B* is the magnetic 
field at is the k component of Bj); Vjj = Vj — is the velocity difference, and Vij^k 

is k component of Vij. The brackets on the left-hand side of the equations mean that the 
physical quantities inside them are evaluated at the position (at the particle i). The term 
ViWij is the gradient of the Kernel function Wij at the position of the particle i, and Hij 
is the artihcial viscosity which allows for appropriate treatment of shock wave dissipation. 
We here adopt the von Neumann-Richtmyer viscosity (e.g., Benz 1990). 


The density and pressure are both calculated from the dehnition of discretness in SPH 
(e.g., Benz 1990; Monaghan 1992), and are expressed, respectively, by the equations: 


N 

(p) = ym,H'y (2d) 

i=i 

N 

(p) = (7 - 1 ) ^ mjUjWij ( 2 e) 

1=1 

As in previous work (GB93; GB94; GB96), the radiative cooling rate (due to collisional 
excitation and recombination) C in equation (Ic), is implicitly calculated using the cooling 
function evaluated by Katz (1989) for a gas of cosmic abundances cooling from from 
T ~ 10® K to T ~ 10^ K. [The cooling is suppressed below T ^ 10"^ K, where the transfer 
of ionizing radiation becomes important, and the assumption of a fully ionized flow breaks 
down (see, e.g., GB93; GBB96)]. 

The choice of a vectorial formalism to write the system of equations above (instead 
of a tensorial form) is to ensure that the magnetic force vector [(V x B) X B], in each 
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point of the system, is strictly perpendicular to the magnetic field vector itself. Tensorial 
implementations in SPMHD, like those proposed by Monaghan (1992), are subject to the 
development of non-physical components of the magnetic force, whose projection on the 
magnetic held vector has non-zero values. Furthermore, these parallel components of the 
magnetic force are proportional to V ■ B (see Brackbill & Barnes 1980). To avoid this, we 
have used the vectorial form and found that those non-physical accelerations are absent in 
our runs [i.e., F ■ B ~ 10“® — 10“’^, where F = ^(V x B) x B]. 

There are several methods which can work well in order to maintain the divergence of 
the magnetic held close to the machine roundoh error in grid-based codes (e.g., Evans & 
Hawley 1988; Schmidt-Voigt 1989). However, the maintenance of V-B = 0 in SPMHD is not 
an easy task. A divergence-cleaning process (e.g., Schmidt-Voigt 1989; Otmianowska-Mazur 
& Chiba 1995), for example, is not applicable yet in the context of the SPMHD, and the 
best way we have found to diminish the ehects of a potential V ■ B 7 ^ 0 was to avoid those 
non-physical accelerations by writing our system of equations in the vector formalism. 
Furthermore, we track the behaviour of V ■ B by evaluating the following quantity at each 
time step and position of the system (e.g, Otmianowska-Mazur & Chiba 1995) 


|V- 

B 

\-h 


|B| 



As in Otmianowska-Mazur & Chiba 1995 (see also Meglicki 1994, 1995 for a detailed 
discussion), we have set u < 10“^ as a limit to the validity of our calculations. In general, 
85% of the particles of the system keep u < 10“^ (with almost all of them distributed 
around the zero value most of the time). For small periods of time, some of them may 
acquire an cu > 10 “^, but the scalar product between the magnetic held and the magnetic 
force remains very small (~ 10 “® — 10 “^, in code units) over the entire evolution of the 
systems simulated here. 
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As additional tests to check the validity of the code modified by the B-£eld 
implementation, we have run Alfven wave tests, similar to those suggested by Clarke (1996). 
In the centre of a rectangular box, a circular pulse of radius r = = 5h (where h is 

the smoothing length) and velocity v = Vinui is set perpendicularly to the magnetic held. 
The latter was assumed to be in the y direction in one case (B = Binit}) and B = Binif{} 
+ k) in the other case. The initial velocity was set to Vinu/vA = ■ Vinit/Binu = 10“^ 

and Vinit/vA = {\/2/2)^/Wp ■ Vinit/Binit = (^2/2)10“^, for each case, respectively (see 
Clarke 1996). Similarly to Clarke’s tests, after several time steps the diffusion of the pulses 
was found to be strictly confined to the direction of propagation of the Alfven waves, as 
depicted in Figs. 17 and 18 of Clarke (1996). 


2.1. Initial and boundary conditions 

The computational domain is represented by a 3-D rectangular box of dimensions 
-17/?, < X < 17/?,, —22Rj < y,z < 22Rj, where Rj is the initial jet radius (/?, is the code 
distance unit). The Cartesian coordinate system has its origin at the center of the box and 
the jet flows through the x-axis, and is continuously injected into the bottom of the box [at 
r = (—17/?j, 0, 0)]. Inside the box, the particles are initially distributed on a cubic lattice. 
An outflow boundary condition is assumed for the boundaries of the box. The particles are 
smoothed out by a spherically symmetric kernel function of width h, and the initial values 
of h were chosen to be OARj and 0.2Rj for the ambient and jet particles, respectively. 

We consider two different initial magnetic held configurations. One of them is an 
initially constant longitudinal B-held permeating both the jet and the ambient medium 
[B = {B^q, 0, 0)]. An observational support for this kind of conhguration is suggested by the 
fact that some protostellar jets appear to be aligned with the main direction of the local 
interstellar magnetic held (e.g., Appenzeller 1989). The other adopted configuration is a 
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force-free helical magnetic field which also extends to the ambient medium with a functional 
dependence given by (see also Todo et ah 1993): 


Br = Q 


(3a) 


B^{r) = Bq 


O.SAdr^ 
.(r+ 0.5d)3 



(3&) 


BJr) = Br 


1 - 


Ar‘^{r + d)] 
(r+ 0.5d)3. 


(3c) 


where r = is the radial distance from the jet axis and the (arbitrary) constants 

A and d are given by 0.99 and 3Rj, respectively. In these equations, Bq is the maximum 
strength of the magnetic held and corresponds to the magnitude of the longitudinal 
component at the jet axis. Figure 1 displays both the longitudinal (B^;) and the toroidal 
{Bip) components as a function of the radial distance (in code units, for = 1 at the jet 
axis). The pitch angle at IRj is ~ 19°. 


The models are parameterized by the dimensionless numbers: i) the density ratio 
between the jet and the ambient medium, y = nj/ua] ii) the ambient Mach number, 

Ma = Vj/ca (where Vj is the jet velocity and is the ambient sound speed); in) the jet to 
the ambient medium pressure ratio at the jet inlet, n = Pj/Pai that we assume to be equal 
to unit; iv) the thermal to the magnetic pressure ratio, at the jet axis f3 = Pth/Pb] and v) 
Qbs = dcooi/Rj, the ratio of the cooling length in the post-shocked gas behind the bow shock 
to the jet radius (see, e.g., GB93). 
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3. The Simulations 

As in previous work (e.g., GB96), based on typical conditions found in protostellar jets, 
we have adopted the following initial values for the parameters: ?] = 3 — 10 (e.g., Morse et 
ah 1992; Raga & Noriega-Crespo 1993), an ambient number density Ua = 200 cm^, Vj = 398 
km s“^ (e.g., Reipurth, Raga & Heathcote 1992), Ma = 24, and Rj = 2 x 10^® cm (e.g., 
Raga 1993). In the MHD simulations, we have assumed an initial (3 = 1 (which corresponds 
to a maximum initial value Bq = 83 fi G). The subsections below present the results of 
the simulations we have performed for both radiatively cooling and adiabatic jets with and 
without magnetic helds, and Table 1 summarizes the values of the input parameters. In 
Table 1, Ma^ = vj/vaj and Mmsj = Vj/ + c|, give the initial Alfven and magnetosonic 
Mach numbers, respectively. The purely hydrodynamical models are labeled with “HD” 
and the MHD models with “ML” and “MH” in the cases of the models with an initial 
longitudinal magnetic field and an initial helical held, respectively. 


3.1. Radiatively cooling jets 

Figure 2 depicts the time evolution of the magnetic held distribution of two 
supermagnetosonic, radiatively cooling jets: ML3r (an MHD model with initial longitudinal 
magnetic held conhguration; top) and MH3r (an MHD model with initial helical magnetic 
held conhguration; bottom). Both models have initial r] = rij/na = 3 and ambient Mach 
number Ma = 24 (see Table 1). We hud that the held lines are amplihed and reoriented 
across the bow shock in both magnetic held conhgurations. The shocked jet material which 
is decelerated at the head and is deposited in the cocoon carries the held lines embedded in 
it. Part of these lines have their polarity reversed as we can see in Fig. 2. 

The cooling parameters behind the bow shock, ~ 8.0, and the jet shock. 
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qjs ~ ri~^qhs ~ 0.3 (see, e.g., GB93) in Fig. 2, imply that, within the head of the jet, the 
ambient shocked gas is almost adiabatic, whereas the shocked jet material is subject to 
rapid radiative cooling. The corresponding density contours and velocity held distribution 
maps of the models above were presented in Figs. 1 and 2 of Paper I, where they were 
also compared with a pure hydrodynamical model (HD3r; see also Table 1). As we have 
stressed in Paper I, the cold dense shell that develops at the jet head due to the cooling of 
the shock-heated jet material in the pure hydrodynamical case, also appears in the MHD 
jets. Similarly, it becomes Rayleigh-Taylor (R-T) unstable (e.g., GB93) and breaks into 
blobs that spill into the cocoon and show a resemblance with the Herbig-Haro objects 
which are observed at the head of protostellar jets. As in the HD case, the density in the 
shell of the MHD jets also undergoes fluctuations with time which are caused by global 
thermal instabilities of the radiative shock (see e.g., Gaetz, Edgar & Ghevalier 1988 and 
GB93 and for details). Fig. 3 shows the time evolution of the shell density {ush) variations 
at the jet axis of the three jets, which have a period of the order of the gas cooling 
time (tcooi ~ 10 years, or tcooi ~ 0.3trf, where td = Rj/ca ~ 38 years corresponds to the 
transverse jet dynamical time). We note that although the HD and the MHD jets attain a 
maximum density of approximately the same magnitude at t/td ~ 1.1, later on the density 
enhancement in the MHD jets is inhibited by the presence of the B-field, particularly in the 
jet with helical held (MH3r). 

Fig. 2 indicates that the cold blobs that develop from shell fragmentation at the head 
of the MHD jets detach from the beam as they are expelled backward to the cocoon. The 
survival of these blobs in the cocoon seems to suffer with the presence of the magnetic 
helds, particularly in the helical case. Gompared to the blobs of the HD jet, their density is 
reduced by a factor ~ 2 in the longitudinal case (ML3r) at the hnal time step {t/td = 1.65), 
and almost vanishes in the helical jet (MH3r). Fig. 4 compares the density and pressure 
prohles across the flow in three different positions along the MHD jets in the head region. 
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The lateral blobs that develop in the cocoon, on both sides of the beam are clearly less 
intense in the helical case than in the longitudinal case. In the helical case, for example, 
we hnd that the toroidal (nonparallel-to-shock) component of the magnetic held, which 
is initially less intense than the longitudinal component (by a factor < 2.5), is strongly 
amplihed by compression in the shocks at the head (by a factor ~ 5). As pointed out in 
Paper I, this amplihcation increases the cooling length behind the jet shock and reduces the 
density growth. As a consequence, the shell tends to stabilize against the R-T instability 
and the clump formation is inhibited. At the end of the evolution, the dense shell with 
clumpy structure which was observed to develop in the head region of the pure HD jet (Fig. 
1 of Paper I) is replaced in the helical case by an elongated plug (Fig. 2) of low density 
material. 

Fig. 2 also indicates the development of some pinching along both MHD jets. In 
the pure HD case, constriction occurs only very close to the jet head where the beam is 
over-conhned by the gas pressure of the cocoon (see Figs. 1 and 2 of Paper I; see also Fig. 
5b below). In the MHD jet with helical held (MH3r), the toroidal component {B^), which 
is amplihed by compression in the shocks at the head, is advected back with the shocked 
material to the cocoon. The associated magnetic pressure (~ B‘^/8n) causes an increase in 
the total pressure of the cocoon relative to the pure HD case, which collimates the beam 
and excites the (fastest growing) small-wavelength pinch modes of the MHD K-H instability. 
These modes over-conhne the beam and drive the approximately equally spaced internal 
shocks seen in the MH3r jet (Fig. 2; bottom). Along the MHD jet with longitudinal held 
(ML3r; Fig. 2; top), the increase in the total conhning pressure of the cocoon also drives 
the development of the MHD K-H instabilities which excite beam pinching and internal 
shocks. Consistent with linear theory for K-H modes, in the supermagnetosonic regime 
considered here, they begin to appear at a distance ~ MmsjRj, which is smaller than in the 
pure HD jet (see Tabel 1). 
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The presence of these oblique internal shocks along the beam of the MHD jets can be 
testhed by the density and pressure prohles across and along the flow depicted in Fig. 5, 
where the pure HD model is also depicted for comparison. We see that the induced internal 
shocks in the MHD cases have a density contrast riis/nj ^ 4 in the ML3r model (middle) 
and Uis/rij 5 in the MH3r model (or Uis/ngh — 0.08 and 0.15, respectively, where Ugh is 
the density at the shell). We also note a close correlation between the pinching zones and 
the appearance of more intense reversed helds in the contact discontinuity between the jet 
and the cocoon. In both magnetic held conhgurations, the original longitudinal components 
are reoriented in the nonparallel shocks at the head and advected back to the cocoon. As a 
consequence, a predominantly toroidal current density distribution (J,^) develops around the 
jet. Such conhguration creates a J x B force (—J^Hn) that constricts the beam triggering 
the MHD pinch formation. 

Fig. 6 depicts the density in the midplane section of the head of three 
supermagnetosonic, radiatively cooling jets with 77 = 10 after they have propagated 
over a distance ~ SSRj: a purely HD model HDlOr (top), an MHD model with initial 
longitudinal magnetic held conhguration, MLlOr (middle), and an MHD model with initial 
helical magnetic held conhguration MHlOr (bottom). Figures 7 and 8 show the time 
evolution of the corresponding velocity and magnetic held distributions. The other initial 
conditions are the same as in Fig. 2. For the time evolution depicted, there is no signihcant 
diherence between the MHD jets with diherent 77 other than those presented by their 
pure HD counterparts. In both HD and MHD cases, the jet with larger 77 [and thus with 
Vbs — Vj /{1 — 77“^/^) — Vj] plows into the ambient medium almost ballistically without 
accumulating much waste of shocked jet gas in the cocoon. The shocked ambient material 
is almost adiabatic {qbs >> 1; see Table 1) and accumulates in the shroud that envelops the 
cocoon/beam, while the shocked jet material cools much faster [qjs << 1 ) and the shock is 
effectively isothermal. The cold shell is thus much thinner in the larger 77 jets and the head 
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resembles a bullet. 

As in the r] = 3 case, the density in the shell of the 77 = 10 jets also undergoes variations 
with time (with a period of the order of the jet radiative cooling time) which are caused by 
the global thermal instabilities of the radiative shocks. Likewise, after reaching a similar 
maximum density amplitude, the MHD jets have their shell density growth inhibited due to 
the decrease of the shock compressibility caused by the presence of the B-field, particularly 
in the jet with helical held (MHlOr). We note, however, that in the HD jet, the shell density 
variations attain a smaller maximum density amplitude than in the MHD cases. This is 
possibly due to the smaller total pressure conhnement that acts on its head. 

In the ?7 = 10 MHD jets (MLlOr and MHlOr), the increase in the total conhning 
pressure in the cocoon (due to the presence of the magnetic held) also over-conhnes the 
beam, relative to the HD jet, and drives some beam pinching. 

Fig. 9a shows the density and pressnre prohles across the how in diherent pinch 
positions along the MHD and HD 77 = 10 jets after they have propagated over a distance 
33Rj (at t/td = 1.65), which can be compared with Fig. 5 (for 77 = 3). As in the 77 = 3 
models, the pinch collimation is larger in the helical MHlOr jet than in the MLlOr jet, 
while the pure hydrodynamic jet, HDlOr, has not developed any internal pinches over the 
time scale depicted. It is interesting to note that in the MHD jets, the fastest-growing 
pinch modes of the K-H instability would be expected to appear only at distances 
MmsjRj 50i?j, which are beyond the computed scales. Thus the early development 
of pinches in these cases, is possibly being triggered mainly by the hydromagnetic 0 -pinch 
effect (e.g., Boyd & Sanderson 1961; Cohn 1983). As in the 77 = 3 case, the jet with 
longitudinal B-£eld also develops a non-axisymmetric helical (kink) mode close to the head 
of the jet which causes some beam twisting. 


As expected, after propagating abont the same distance, the pinches that develop in 
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the larger rj case are less numerous, since the total amount of conhning shocked material 
that deposits into the cocoon is comparatively smaller (see Figs. 5 and 9 for a comparison). 
Nonetheless, the density contrasts, riis/rij, attained in the pinch regions of the MHD jet 
with longitudinal held, MLlOr {riis/rij ~ 3 — 4.5, or riis/nsh ~ 0.05 — 0.07) are approximately 
the same as those in the smaller rj jet (ML3r). The situation is a little more complex for 
the MHD jets with helical held. Over the whole jet evolution, some few pinches are found 
to be stronger in the larger rj jet (MHlOr) (Fig. 8), but in general the density contrasts are 
approximately the same in both cases {riis/rij ^ 5, or riis/nsh ^ 0.13). 


3.2. Adiabatic jets 

Fig. 10 depicts the density in the midplane section of the head of three 
supermagnetosonic, adiabatic jets with rj = 3: a purely HD jet (HD3a, top), an 
MHD jet with initial longitudinal B-held, (ML3a, middle), and an MHD jet with initial 
helical B-held (MH3a, bottom). Figures 11 and 12 show the time evolution of the 
corresponding velocity and magnetic held distributions. The initial conditions are the 
same as in Fig. 2 (see Table 1). Previous numerical analyses comparing pure HD jets 
with and without radiative cooling (e.g., BFK, GB93) have shown that the presence of 
radiative cooling tends to reduce the strength and number of internal shocks excited by K-H 
instability. Consistently, the HD3a jet in Figs. 10-12 reveals the appearance of a pinching 
zone (at x ~ IRj) which is absent in its radiatively cooling counterpart (HD3r, Fig. 2; 
see also Figs. 1 and 2 of Paper I). Besides, the beam constriction which appears close to 
the head is larger in the adiabatic jet (HD3a). These results are also compatible with 
the predictions of the linear stability theory (Hardee & Stone 1997), when applied in the 
context of the cooling function employed in this work (see §1 ). Similarly, the pinches which 
develop in the MHD adiabatic jets (Figs. 10-12) are generally more intense and numerous 
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than in the radiatively cooling counterparts of Fig. 2. Their larger strength can be testihed 
by direct comparison of Figs. 13 and 5 which show the transverse prohles of some pinches 
along the jets. In particular, for the ML3a jet, we hnd a chain of 8 evolving pinches, at 
t/td = 1.65, against 5 in the radiatively cooling jet (ML3r), and the corresponding densities 
nis/rij, are about 2 to 3 times greater than in the ML3r jet. The helical adiabatic jet MH3a 
displays pinches that are 40% denser than those found in the radiatively cooling jet (MH3r) 
at t/td = 1.65. 

We hnd similar results when comparing radiatively cooling and adiabatic jets with 
larger rj. Figure 14 shows the midplane density (left) and velocity held distribution (right) 
for adiabatic iets with 77 = 10 after they have propagated ~ 33Rj, with the same initial 
conditions as the cooling jets of hgures 6 , 7, and 8 . The corresponding magnetic held 
distribution is presented in Fig. 15. As before, the 77 = 10 adiabatic jet with longitudinal 
B-held (MLlOa) develops stronger pinches than its cooling counterpart (MLlOr model) by 
a factor ^ 3. The 77 = 10, adiabatic helical jet (MHlOa), on the other hand, has pinches 
with densities of the same order of magnitude as those in the cooling jet (MHlOr). Also, 
the pinches are found to be stronger in the adiabatic jets with smaller 77 , in both magnetic 
held conhgurations. 


4. Conclusions and Discussion 

We have investigated here the ehects of magnetic helds on the structure of evolving 
overdense, radiatively cooling, supermagnetosonic jets with the help of 3-D SPMHD 
simulations and compared with purely hydrodynamical and adiabatic calculations. Two 
initial magnetic held conhgurations (with magnitude in approximate equipartition with 
the gas) have been examined: a longitudinal, and a helical held permeating both the 
jet and the ambient medium. Calculated for a set of parameters which are particularly 
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appropriate to protostellar jets (with density ratios between the jet and the ambient medium 
f] ^ 3 — 10, and ambient Mach number ^ 24), our results indicate that magnetic helds 
have important effects on the dynamics of radiatively cooling jets. Both magnetic held 
geometries are able to improve jet collimation relative to the pure hydrodynamical (HD) 
jets, but this effect is larger in the helical held case. 

As we have stressed in Paper I, the cold dense shell that develops at the jet head 
due to the cooling of the shock-heated jet material in the HD cases, also appears in the 
MHD jets. Likewise, it becomes Rayleigh-Taylor (R-T) unstable and breaks into clumps 
which are more visible in the smaller rj jets where the developed shell is thicker. Also as in 
the HD case, the shell of the MHD jets undergoes density variations with time which are 
caused by global thermal instabilities of the radiative shocks. However the amplihcation 
and re-orientation of the nonparallel components of the magnetic helds by the shocks at 
the head, particularly in the helical held geometry, reduces the postshock compressibility 
and increases the postshock cooling length. This tends to stabilize the shell against both 
the R-T and the thermal instabilities. As a consequence, the clumps that are observed 
to develop by fragmentation of the shell in the HD jets are depleted in the helical held 
geometry. The jet immersed in the longitudinal held, on the other hand, still retains the 
clumps although they have their densities decreased relative to the the HD counterparts. 
The fact that the clumpy shell structure resembles the knotty pattern of the Herbig-Haro 
objects which are commonly observed at the head of protostellar jets (e.g., Herbig & Jones 
1981; Brugel et ah 1985; Reipurth 1989; Heathcote et al. 1996) suggests that a longitudinal 
magnetic held geometry would be more likely in the outer regions of these jets than a helical 
held geometry (see also Paper I). 

Over the computed time and length scales, internal oblique shocks along the beam are 
not found to develop in the HD systems examined here. On the other hand, in their MHD 
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counterparts with both magnetic held conhgurations, the conhning total pressure of the 
cocoon excites (the fastest growing) low-amplitude MHD K-H rehection pinching modes 
which drive a chain of approximately equally spaced internal shocks along the beam, but 
these shocks are found to be slightly stronger in the helical held case (by a factor 20%). 
Also, as expected, the internal shocks tend to appear in a larger number in the smaller 
rj jets (due to the larger amount of conhning shocked material which is deposited in the 
cocoon), although their densities are approximately of the same magnitude as those in the 
larger rj jets. A non-axisymmetric helical mode is also excited close to the head of the 
radiatively cooling MHD jets with longitudinal held, which causes some beam wiggling. 

The number and strength of the internal shocks excited in the MHD adiabatic jets is 
larger than in the radiatively cooling counterparts (by a factor ^ 2 in both number and 
strength). This result is compatible with the linear stability theory (Hardee & Stone 1997) 
when applied in the context of the cooling function employed in the present work, and is 
also compatible with previous numerical work of HD jets which has shown that the presence 
of radiative cooling tends to reduce the strength and number of internal shocks along the 
jet (e.g., BFK; GB93). Also, the pinches are found to be more numerous in the adiabatic 
jets with smaller rj, in both magnetic held conhgurations. 

The internal shocks are found to propagate downstream with velocities close to that 
of the jet head {vbs ~ 250 km s“^). The mean distance between them (^i 2 — 47?^) is in 
agreement with the observed knots in the jets. However, the weakness of the shocks in 
the radiatively cooling jets {jiis/nj ~ 3 — 5) makes it doubtful that they could produce 
by themselves the bright knots observed in protostellar jets. Probably, other mechanisms, 
like intermittency in the jet injection velocity, play a more relevant role in knot formation 
in those jets (see e.g., Raga et al. 1990; de Gouveia Dal Pino and Benz 1994; Stone and 
Norman 1993a). 
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We should note that in the recent numerical study of magnetized radiatively cooling 
jets by Frank et al. (1998), they have found that toroidal helds (also in approximate 
equipartition with the gas) may excite the development of strong pinches along the beam. 
This apparent contradiction between their analysis and ours is possibly due to differences 
in the assumed initial conditions. Frank et al. have considered a slower and much lighter 
jet (with a jet Mach number Mj = Vjjcg. ~ 10 and ri = 1.5) than the cases examined here 
(Mj = Vjjcsj ^42 — 76 and rj = 3 — 10). Thus, consistently with our results above, their 
smaller r] jet should be expected to produce a larger amount of more intense pinches. This 
result has also been conhrmed by numerical simulations (not presented here) that we have 
performed of ?7 = 1 jets, which have produced a larger amount of pinches with slightly 
larger densities than the larger rj jets studied above. 

Finally, we should make some remarks on the late evolution of the radiatively cooling 
magnetized jets. In the magnetic held maps of Figs. 2, and 8 above, we have detected the 
development of (sometimes strong) magnetic held reversals at the contact discontinuity 
between the jet and the cocoon with intensities up to 5 times their initial magnitude. 

As stressed in §3, held reversals occur in both investigated magnetic held conhgurations, 
because the held lines are amplihed by compression in the nonparallel shocks at the jet 
head, and are enforced to how backward with the shocked plasma into the cocoon. In this 
process, the lines are reoriented and sometimes have their polarization reversed. Beyond 
the integrated length and time scales depicted in the hgures above, however, the increasing 
strength of the reversed helds due to shear at the contact discontinuity (see e.g., eq. 

2 c), may lead to the development of strong pinching regions which ultimately may cause 
jet disruption, particularly in the cases with longitudinal held. As an example. Fig. 16 
depicts the time evolution of the velocity (left) and magnetic held (right) distributions of 
a disruption zone which occured in the late evolution of the ML3r jet of Fig. 2. We can 
clearly distinguish two regions with reversed B-helds whose strength increases with time. 



which are correlated with developing pinches. The inner constriction becomes so strong 
that it finally causes the disruption of the beam. Although shear and compression are 
expected to enhance B, this amplification is possibly partially due to numerical effects. 

In fact, along a contact discontinuity with such a magnetic field topology with oppositely 
directed field lines pressed together, magnetic diffusion and reconnection may have an 
important role and lead to intense magnetic energy release. Of course, under the ideal-MHD 
approach considered here, these dissipation effects of the magnetic field have not been 
appropriately considered thus leading to possibly anomalous amplification of the reversed 
components in the late stages of the evolution of some of the jets. Further, we should 
note also that no jet disruption was detected in the majority of the adiabatic cases. This 
fact is consistent with recent numerical studies of reconnection processes in 2-D current 
sheets (Oreshina & Somov 1998), which have shown that reconnection rates are smaller in 
adiabatic than in radiatively cooling plasmas. Although the transposition of these results 
to the more complex flow geometry we have investigated here is not straightforward, they 
seem to suggest that the appropriate consideration of magnetic field dissipation in our 
models will possibly decrease and even suppress the disruptive effects of the magnetic 
fields found in some of the radiatively cooling cases examined here in their late evolution. 
The transformation of magnetic energy into thermal energy of the gas will probably have 
important effects on the structure and emission mechanisms of the beam (see e.g., Malagoli, 
Bodo & Rosner 1996; Min 1997a, b; Jones et al. 1997) and also in the process of turbulent 
mixing of the jet and cocoon material. Such finite magnetic resistivity effects will be 
addressed in a forthcoming paper. Besides, the potential signatures that magnetic fields 
may leave on the morphology of radiative cooling jets, especially behind the shocks, provide 
important constraints that can be used in future observational tests to distinguish among 
different candidate mechanisms for emission line production, jet collimation, and turbulent 
entrainment at the contact discontinuity between the jet and the cocoon. 
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TABLES 

Table 1: The models and its initial physical parameters. 


Run 

Ma 

Ma- 

Mms, 

h 

/9 

Qjs 

Qbs 

HD3a 

24 

— 

41.6 

3 

oo 

— 

— 

HDlOa 

24 

— 

75.9 

10 

oo 

— 

— 

ML3a 

24 

38.1 

28.0 

3 

1 

— 

— 

MH3a 

24 

38.1 

28.0 

3 

1 

— 

— 

MLlOa 

24 

70.6 

51.4 

10 

1 

— 

— 

MHlOa 

24 

70.6 

51.4 

10 

1 

— 

— 

HD3r 

24 

— 

41.6 

3 

oo 

0.29 

8.03 

HDlOr 

24 

— 

75.9 

10 

oo 

0.02 

16.56 

ML3r 

24 

38.1 

28.0 

3 

1 

0.29 

8.03 

MH3r 

24 

38.1 

28.0 

3 

1 

0.29 

8.03 

MLlOr 

24 

70.6 

51.4 

10 

1 

0.02 

16.56 


MHlOr 24 70.6 51.4 10 1 0.02 16.56 
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FIGURE CAPTIONS 

Figure 1: Longitudinal {Bx'-, left) and toroidal right) magnetic field components (see 
equations 3 in the text) as a function of the radial distance [r = + 2 ;^)^/^], for /? = 1 at 

the jet axis. The coordinates are in code units (for the magnetic helds 1 c.u. = 21.5/i G). 

Figure 2: Midplane magnetic field distribution evolution of the 77 = 3 radiatively cooling 
jets ML3r (top) and MH3r (bottom). The initial conditions are r] = rij/ua = S, Ua = 200 
cm“^, Ma = 24, Vj ~ 398 km s“^, /? = Sirp/B'^ ~ 1, ~ 8 and qjg ^ 0.3. The times 

and the jet head positions are: t/td = 1.40 and x ~ 25i?j (left); t/td = 1.60 and x ~ 29i?j 
(middle); t/td = 1.65 and x ~ 30i?j (right). Note the reorientation (and amplihcation), 
in both models, of the magnetic helds that are carried with shocked jet material into the 
cocoon. 

Figure 3: Density of the shell at the jet axis as a function of the time for the three 
radiatively cooling jets presented in Fig. 2: HD3r (solid line); ML3r (dotted line) and MH3r 
(dashed line). 

Figure 4: Density (solid line) and pressure (dashed line) across the how at diherent 
positions along the how in the head region of the ML3r (top) and MH3r (bottom) jets. The 
positions are in units of Rj. The central peak corresponds to the beam region, while the 
secondary peaks on both sides of the beam are the blobs which are much denser in the ML3r 
model (top) than in the MH3r model (bottom). These prohles are taken at t/td = 1.65, and 
the density and pressure scales can be calibrated using the marker in the top region of the 
plots. 
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Figure 5: (a) Density (solid line) and pressure (dashed line) profiles across the HD3r 
(top), ML3r (middle), and MH3r (bottom) jets taken in three different pinching positions 
along the flow (at t/td = 1.65). All the prohles have been scaled as in Fig. 4. (b) The 
corresponding axial density along the jet axis {y = z = Q) showing the channel of internal 
shocks. The position of the three pinches depicted in (a) are labeled with arrows. Note that 
the origin of the coordinates along the jet axis has been shifted from —QRj to 0. 

Figure 6: Gray-scale representation of the midplane density of radiative cooling jets with 
?7 = 10: a hydrodynamical jet (top; model HDlOr), an MHD jet with initial longitudinal 
magnetic held distribution (middle; model MLlOr) and an MHD jet with initial helical 
magnetic held distribution (bottom; model MHlOr), at a time t/td = 1.65. The initial 
conditions are y = Uj/ua = 10, Ua = 200 cm“^, Ma = 24, Vj ~ 398 km s“^, q^s — 16 
and Qjs — 0.02. In the MHD cases, the initial /5 ~ 1. The gray scale (from minimum 
to maximum) is given by black, light-gray, white and dark-gray. The maximum density 
reached by the shell at the head of the jets is Ugh/na — 750 (top), Ugh/na — 550 (middle) 
and Ugh/Ua — 450 (bottom). (This and the remaining gray-scale plots below were built with 
the help of the GREY subroutine developed by Marinho & Lepine 1998). 

Figure 7: Midplane velocity held distribution evolution of the jets of Fig. 6: HDlOr 
model (top), MLlOr model (middle) and MHlOr model (bottom). The times and the jet 
head positions are t/td = 1.40 and x ~ 27i?j (left); t/td = 1.60 and x ~ ?>lRj (middle); 
t/td = 1.65 and x ~ 33i?j (right). 

Figure 8 : The same as in Fig. 7 for the magnetic held distribution of the MLlOr model 
(top) and MHlOr model (bottom). 
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Figure 9: The same as in Fig. 5 for the rj = 10, radiatively cooling jets: HDlOr (top), 
MLlOr (middle), and MHlOr jets (bottom), at t/td = 1.65. 

Figure 10: Gray-scale representation of the midplane density of adiabatic jets with t] = 3: 
a hydrodynamical jet (top; model HD3a), an MHD jet with initial longitudinal magnetic 
held distribution (middle; model ML3a) and an MHD jet with initial helical magnetic 
held distribution (bottom; model MH3a), at a time t/td = 1-65. The initial conditions are 
the same as in Fig. 6. The gray scale (from minimum to maximum) is given by black, 
light-gray, white and dark-gray. The maximum density reached at the head of the jets is 
nsh/ria — 37 (top), Ugh/na — 39 (middle) and Ugh/na — 42 (bottom). 

Figure 11: Midplane velocity held distribution evolution of the jets of Fig. 10: HD3a 
model (top), ML3a model (middle), and MH3a model (bottom). The times and the jet 
head positions are t/td = 1-40 and x ~ ‘^7Rj (left); t/td = 1-60 and x ~ SlRj (middle); 
t/td = 1-65 and x ~ 33Rj (right). 

Figure 12: Midplane magnetic held distribution of the adiabatic MHD models presented 
in Fig. 11: ML3a (top) and MH3a (bottom). 

Figure 13: (a) Density (solid line) and pressure (dashed line) prohles across the HD3a 
(top), ML3a (middle), and MH3r (bottom) jets taken in three diherent pinching positions 
along the how (at t/td = 1-65). The prohles have been scaled as in Figs. 4 and 5. (b) The 
corresponding axial density along the jet axis {y = z = R), showing the channel of internal 
shocks. The position of the three pinches depicted in (a) are labeled with arrows. 


Figure 14: Gray-scale representation of the midplane density (left) and velocity held 
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distributions (right) of adiabatic jets with rj = 10: a hydrodynamical jet (top; model 
HD 10a), an MHD jet with initial longitudinal magnetic held distribution (middle; model 
MLlOa), and an MHD jet with initial helical magnetic held distribution (bottom; model 
MHlOa), at t/td = 1.65. The initial conditions are the same as in Fig. 6. The gray scale 
(from minimum to maximum) is given by black, light-gray, white and dark-gray. The 
maximum density reached by the shell at the head of the jets is Ugh/na — 182 (top), 
nsh/na — 170 (middle) and Ugh/na — 140 (bottom). 

Figure 15: Midplane magnetic held distribution for the jets of Fig. 14: MLlOa (top), and 
MHlOa (bottom) at t/td = 1.65. 

Figure 16: Time evolution of the velocity held (left), and the magnetic held distribution 
(right) of the ML3r disrupting region. The times shown, from top to bottom, are: 
t/td = 1.40, 1.45, 1.49, 1.52 and 1.55. See the text for discussion. 
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